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Abstract. The transverse spin fluctuations are introduced to the density functional theory for supercon- 
ductors (SCDFT). Paramagnons are treated within the random phase approximation and assumed to be 
the same for the normal and superconducting state. The effect of spin fluctuations on T c is studied for a 
few simple metals at ambient pressure and niobium at several pressures up to 80 GPa. 
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1 Introduction 



Since the discovery of superconductivity many theories 
have been born to explain this phenomenon and calcu- 
late observables. First papers about the role of spin fluc- 
tuations by Doniach and Engelsberg pQ and Izuyama et 
al. were published in sixties. Till today, fluctuations 
have been introduced to many-body and phenomenologi- 
cal models and a very popular semiempirical theory pro- 
posed by Eliashberg 0. 

The goal of this work is to include the spin fluctu- 
ations into the density functional theory for supercon- 
ductors which, in principle, anables to calculate all ma- 
terial properties, also in the superconducting state, from 
first principles. The framework of the SCDFT was set up 
by Oliveira, Gross and Kohn 0] in 1988. Recently, the 
SCDFT gap equation has been solved numerically for sim- 
ple metals [HUE] and MgB 2 

As for the critical temperatures, it is known for a long 
time, that spin fluctuations decrease considerably T c of 
some superconductors |HKj ■ In our previous work for nio- 
bium under pressure |1(J| . we solved the gap equation of 
the Eliashberg theory with and without spin fluctua- 
tions and the SCDFT gap equation only with the Coulomb 
and phonon interactions. We found that the effect of para- 
magnons decreased T c obtained from the Eliashberg the- 
ory by 3-4 K, however, an approximate treatment of the 
Coulomb interactions by a simple constant, fx*, led to a 
large disagreement of the theoretical results with the ex- 
perimetal data (TT]. In contrast to the Eliashberg theory, 
the SCDFT scheme is parameter free, but the critical tem- 
perature calculated without spin fluctuations for Nb at 
ambient pressure JU| was about 3.7 K higher than the 
experimental T c . 



In this work, we follow the derivations of the SCDFT 
gap equation given in a number of PhD theses 1 |12U13lll4j , 
and we include the spin fluctuations. The paramagnon 
spectral function is calculated within the random phase 
approximation (RPA) with the assumption of the homo- 
geneous electron gas, similarly to the work by Berk and 
Schrieffer |S] done for the Eliashberg theory. We solve the 
obtained gap equation for a few simple metals and update 
our previous results for niobium under pressure. 

In the following Sections, we introduce the SCDFT gap 
equation and the construction of the exchange-correlation 
functional, F xc , by collecting the most important building 
blocks of the theory given by its authors 0j and first de- 
velopers [t21|13II14II15II16| . These Sections are: II. SCDFT 
gap equation, III. Exchange-correlation functional, and 

IV. Coulomb interaction and phonons in F xc . Above Sec- 
tions are written using the notation according to Parks 
[T7irr%| and Vonsovsky \nj\. This notation is in some points, 
such as Nambu Green's function and the selfenergy, dif- 
ferent than the notation previously used for the SCDFT 
[T^lll3lll4j . We introduce the spin fluctuations in Sections: 

V. Paramagnons in F xc and VI. Gap equation with para- 
magnons and implementation details. We report obtained 
critical temperatures in Section VII, and we summarize in 
Section VIII. 



2 SCDFT gap equation 

In this Section, we wish to guide the reader, step by step, 
to the gap equation which will be solved at the end of 
this work to calculate the critical temperatures. We start 
by bringing the fundaments of the SCDFT £|] and the 
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main approximations, such as the decoupling of band en- 
ergies and the superconducting gap and a linearization of 
the gap equation close to T c , which were assumed for nu- 
merical convenience |12H13U14II15| . We believe that these 
approximations do not cause any significant difference in 
the calculated critical temperatures. 

Turning to details of the SCDFT, in order to obtain 
the gap equation one needs to follow the points below: 

1. The grand- canonical Hamiltonian for a superconductor 
reads 



Hv,A = J2 



+ v(r) - \l 



Vv(r) 



+ i f d 3 r d 3 r' &(r)#( r ') —L_ ^(r')^(r) 



d 3 n d 3 r[ d 3 r 2 d 3 r' 2 ^{(^{(n) 
x w{r' 1 ,Y U Y 2 ,v' 2 ) ^ T (r2)^4.(r^) 

d 3 r d 3 r' A* (r, r') ^ (r)tp i (r') + H.c. 



(1) 



where v(r) and Z\(r,r') are an external potential and 
an anomalous pair potential respectively. The pair- 
ing interaction w in the particular BCS case satisfies 
w(r[, ri, r 2 , t' 2 ) = w(r[ — ri, r 2 — r 2 ). The normal and 
anomalous densities, n(r) and x( r J r ')i are defined as 



n(r)= ^{^(^(p)), 
X (r,r') = (W^l(r')). 



(2) 
(3) 



2. The Hohenberg-Kohn theorem for superconductors says 
that, at each temperature 6 = 1/(3, the normal and 
anomalous densities, n(r) and x(r, r'), determine uniquely 
the density operator p — e~^ Hv - A /Tre~^ Hv - A which 
minimizes the thermodynamic potential, f2 v a [p] , given 

by 

f2 ViA [p] = Tr{pH v<A + e plnp}. (4) 

3. Furthermore, the thermodynamic potential can be ex- 
pressed in terms of the densities and potentials by in- 
volving a universal functional of the densities, F[n, x] , ^ 
as follows 



&v,A[n,x\ = F[n, x] + J d 3 r v(r)n(r) 

- J d 3 rd 3 r'[A*(r,r') X {r,r') + H.c.].(5) 

4. The universal functional contains the exchange- correlation 
(xc) free- energy functional, F xc [n, x], as below g 

F[n,x]=T s [n,x]-eS s [n, X ]-^N 

+ 1 fd 3 r d 3 r' n(r) " (r,) 
+ 2 J drdr | r _ r /| 

- J d 3 n d 3 r[ d 3 r 2 d 3 r' 2 X*(n,ri) 

x u;(ri,ri,r 2 ,r 2 )x*(r2,r 2 ) + F xc [n, X ], (6) 



where T s [n, x] an d S s [n, x] are the kinetic energy and 
the entropy of a noninteracting system with the nonin- 
teracting potentials, v s and A s , such that the densities 
n and x are equal to those of the noninteracting sys- 
tem. In the above formula, p is the chemical potential. 
The noninteracting grand- canonical Hamiltonian can 
be written in terms of the noninteracting densities and 
potentials as 



H„ = 



+ v s (r) - p 



d s r d 3 r' A* (r, r')V> T (r)V>j. (r') + H.c. 



■(7) 



The diagonalization of the noninteracting Hamiltonian, 
H s , using the Bogoliubov transformation leads to the 
Kohn-Sham-Bogoliubov-de Gennes (KS-BdG) equations 

V '~ 1 Ui(v) + [ d 3 r' A s {v,v') Vl {v') 



2 

V 2 



+ v s (v) - p 
+ v s (r) - p 



= Ei Ui(r), (8) 
»-(r) + J d 3 r' zl*(r,r'K(r') 

= Ei Vi (r), (9) 

with Vi(r) and Ui(r) being the pair creation and ani- 
hilation amplitudes respectively. 

The noninteracting potentials, v s and A s , consist of 
the external potentials, vq and Aq, and Hartree po- 
tentials, and the exchange-correlation potentials, v xc 
and A xc , as follows 



v s [n, X }(r)=vo(r)+ I d 3 r' ^1 



+ v xc [n, x](r), 
A s [n,x]{r,r') = A {r,r') + J 
+ A xc [n,x]{v,v') 



(10) 



d 3 r'^ll 



(11) 

The external pairing potential has been introduced in 
order to break the symmetry, thus, in the calculations 
A)(r,r') — > in Eq. dTTJ) . 

The exchange- correlation potentials, v xc and A xc , are 
defined as the derivatives of the xc functional, F xc [n, x], 
with respect to the densities, n and x> correspondingly 
as below 

SF xc [n,x] 

v X c[n,x}(r) = , , - , (12) 



4jrc[ra,x](r,r') = - 



8n(r) 
SE xc [n,x] 



(13) 



£x*(r,r') ' 

The densities, n and x are defined as functions of the 
amplitudes Uj(r) and Uj(r) as 

n(r) = 2^ [ + \vi(r)\ 2 (l ~ fai) }, (14) 



C(r,r') = E [ v*(r') Ui (r)(l - fo) v* (r)u t (r')f^ 



(15) 
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with the Fermi distribution function fpj = l+exp((3Ei). 

At this point, one could guess the densities, n and x, 
and find the potentials, v xc and A XCy and solve the KS- 
BdG equations, and find new densities etc. Further for 
practical reasons, as we already mentioned at the begin of 
this Section, one can make two approximations which we 
will discuss now. 

10. The energy scales for the electronic energies and the 
superconducting energy gap differ by orders of magni- 
tude. Therefore, the KS-BdG equations can be decou- 
pled into the Kohn-Sham equation and the gap equa- 
tion. This approximation was introduced to the SCDFT 
in Ref. Q2]. 

It holds within the decoupling approximation that: 

(a) the amplitudes iti(r) and Vi(r) can be written in a 
form 

itj(r) w m <pi(r) ; (r) « Vi <fi (r) , (16) 

(b) the eigenvalues in Eqs. JHJ and © are defined by 

Ei = ±y/& + \A (17) 

where & = £j — fi 

(c) the coefficients m and Vi are given by 



The Eq. (jJEJ) stems from including Eqs. and (|2"3")l 

into Eq. I|ll|l , and using the potential given by formula 
(fTT|> in Eq. jST) , 

12. In vicinity of T c , the gap function is vanishing, there- 
fore, it can be linearized in Aj. 

The above twelve steps lead to the gap equation which 
can be expressed in the form 

A = -iE^-M^|^A-, (26) 



M H xc,ij [/•*] 



5A 



Hxc.i 



(27) 



±1***)* ^1 + | 



1 



(18) 
(19) 



(20) 



1 

Vi " 71 v ■ ^ 

and the phase factor (f>i is defined by 

(d) the matrix elements Aj are defined as 

A i =fd 3 r J d 3 r' ^*(r)A(r,r> 4 (r'), (21) 

(e) and the normal and anomalous densities read re- 
spectively 



n(r)=£(l 



x(r,r')4E 



-J|)tanh(^) r ,(r.h'. (22) 
|i t anh^)^(r)^(r'). (23) 



11. The decoupling of the two energy scales yields a trans- 
formation of the KS-BdG equations into the ordinary 
Kohn-Sham equation 



V 2 



+ v s [n,x}(r) - n 



<Pi{r) = e t ifiiv), (24) 



and the gap equation 

Ai = A Hx ci[n, Ai]. 



where Ah xc ,i is defined by Eq. (|13fl . 
In other way, Eq. (|27|l can be written as 



A = -ZtMAt -l^KvW^&Ai. (28) 

Kij and Z t are the functionals only of the chemical poten- 
tial in the case when the gap equation is linearized. The 
above gap equation will be solved later in this work. The 
explicit form of the kernel and the norm Zi will be 
given in Sectional 

Since the gap function (|25l) contains the exchange- 
correlation part defined by Eq. (|13|) . we will focus on the 
construction of the exchange-correlation free-energy func- 
tional, F XCl in the following Section. 



3 Exchange-correlation functional, F xc [n,x\ 

The derivation of the exchange-correlation energy F xc , by 
making use of the perturbative expansion of the thermo- 
dynamic potential, was given in Ref. |12j . For the purpose 
of inclusion the spin interactions, we will briefly draw a 
skeleton of this derivation here. 

First, one can notice from Eqs. JHJ) and © that 



d 3 r [v H (r) + u xc (r)]n(r) 



- J d 3 rd 3 r' [A* xc (r,r')x(r,r')+A xc (r,r') X *(r,r')} 

(29) 



- / d 3 r d 3 r' n ( r ) n ( r ') , 
2.7 r-r' 



Then, we take the coupling constant integration for- 
mula which reads 



1 dX 
o ~' 



(30) 



(25) 



where A is the coupling constant, and the perturbation 
Hamiltonian Hi satisfies H = H s + XHi with the inter- 
acting and noninteracting Hamiltonians, H and H s , re- 
spectively. The Hamiltonian Hi contains the difference 
between the exact Coulomb interaction and the exchange- 
correlation potentials, the electron-phonon interaction, the 
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electron- paramagnon interaction etc. _ ^ / ( j3 r \ VH ( r \ _|_ v xc (r)] n x (r) 

The average in Eq. (|30|) has to be taken with the density J 

operator p x = e~ 0Hx /Z\. [331* 1 \ i\ 

Before an explicit evaluation of the coupling constant + 2X d r d r A xc (r,r )\ (r, r ) > , (41) 
integration formula (|3U[1 , we write here a definition of the 

Nambu Green's function which we can plug into the Eq. for the exchange- 



es (rr, r'r') = 



correlation functional, F xc [n,x]- 



Gee' (rr, r r ) F a -„i (rr, r r ) \ As (or the first- order selfenergy, £ a &' , ioi the nonm&g- 



F-ee'( rT ' r ' T ') — G_ CT /_ ff (r r , rr) y ' netic systems with the potential u(rr, r'r'), this energy is 

defined as 

which is a 2 x 2-matrix of the normal and anomalous single 

particle Green's functions, G aa > and F aa > , given respec- S(tt, r'r') = —v (rr, r'r') TaGlrr, t't')ts, (42) 

tively by 

G CTCT ,(rr,rV) = -(Ti(rr)^(rV')), (32) and T3 iS ° ne ° f the Pauli matriC6S: 
F CTCT ,(rr,r'r') = -(T^(rr)Vv (r'r')>, (33) ^ = ( A ^JOi 

Fl,{vr,v'r') = -{ni(rr)^Av'r')). (34) V 10 /' V"* 

_ 1 A (10 

The detailed derivation of {H x ) is given in Refs. [121101 T3 ~ ^ -1 ) ' T ° ^01 

I14| . This derivation starts from the equations of motion for 
the field operator, Vv, and for the noninteracting Green's For the magnetic systems, the matrix r 3 in the each vertex 
function, G^ CT , , which are as follows of Feynman diagrams for the selfenergy with the Coulomb 

and phonon interactions has to be replaced with the ma- 



^-MrT)=e^[H,i, a (r)]e-^, (35) 



trix r T 3 . 



t G s aa , (rr, r'r') = — S aa i 5(r — r') <5(r — r'), (36) In this Section, we sketched main steps to be done 

for finding a general form of the F xc [n, x] functional for a 
with the Kohn-Sham Hamiltonian for the normal state, superconductor. The final formula involves the selfenergy 
h s , and the operator £ given respectively by which will be evaluated in detail for the Coulomb and 

electron-phonon interactions in the next Section and for 

hJr) = -— + v s (r) - \i (37) tne P arama g nons in Section 

d 



(38) 

A*(r) £p -h s (r)J 4 Coulomb and electron-phonon interactions 

The operator A s (r) is defined as ' xc iX 

/ 3 , - , , The derivation of F xc for the Coulomb and phonon inter- 

d r Z\ s (r,r )/(r ). (39) actions is given in detail in Refs. |51ll3| . Here, we report 

this derivation starting with the interactions in the self- 

In order to complete the derivation, one also needs to make energy (in Eq. 142|) ') defined by 
use of the Dyson's equation 

G^(rr,r'r') = G*, CT (rr,rV) v° l (r,r>) = j-L^, (43) 

+ E/ d ^ d ^ J dr i dr 2 ^(rr.rm) v ph (rr,r'r') = F Aq (r)£ Aq (r - r')l/ Aq (r'), (44) 

ee' 

x I7(riri Y2T2) G a e'(j2T2 r'r') (40) where V\q is the electron-phonon interaction vertex and 

D\q is the phonon Green's function defined as 

with S being the self-energy. 

£> Aq (r,r') = (f^ q (r)^ q (r')), (45) 

The above building blocks make us to arrive, after 

some algebra, at the relation ■ ,i , , u t , ,t /> \ 1 ■ a 1 

with 3> Aq = 6 Aiq + o' x , and b' x (o A , q ) being the phonon 

1 fi j\ f r creation (anihilation) operators. 



Let us have a look now at the expression (|29|l for F xc 

/ f _ A 1 1 -\ 11 + an( I the definitions of the Nambu Green's function and 

di~ [£ee> ( rT > r T )G a i a (r t ,rr )] u selfenergy given by Eqs. (jSU a nd JHJ| respectively. The 
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(l,l)-element of the (i7G)-matrix, present in the formula 
l|41(l and enterring Eq. (|29|l . is proportional to 



It 



(46) 



and the corresponding terms with the opposite spins. The 
above terms appear for both the Coulomb and electron- 
phonon interactions, and later will lead to the opposite 
signums in the kernel K%j and the norm Z% of the gap equa- 
tion. Just mentioned difference in signum, in the first order 
terms of the total energy with the normal and anomalous 
Green's functions, stems from the factor of (-1) which one 
has to associate with the each loop of anomalous Green's 
functions. 

In order to evaluate further F xc , we bring here the 
explicite expressions for the noninteracting propagators. 
The formulas given below were derived from the defini- 
tions I|32I34[) assumming the decoupling approximation, 
i.e. Eqs. (|16fl : the Kohn-Sham orbitals y>k(i") were chosen 
to those of a homogeneous gas (w n are the odd Matsubara 
frequencies) 



G s aa ,(k,w n ) = S^ 



_iuj n -E k iui n + E k 
FJ CT ,(k;w;„) = S a - a > sgn(a') 

* ( 1 1 
x u kv k - — =- - 

\ lU) n + E k lUJ n - E k 

F*l,(k;w n ) = 5 a - a > sgn(a) 

* ( 1 1 

\ lLU n + E k lUJn - E k 



(47) 



(48) 



(49) 



Now, we will combine Eqs. lj29~jl and for the F xc 

and Q — C2 S respectively, with a definition of the Nambu 
Green's function, Eq. and an expression for the self- 
energy, Eq. H42|l. As for the noniteracting Green's func- 
tions, we use those obtained within the decoupling approx- 
imation, i.e. 14714911 . This way, one arrives to formulas for 
the xc energy, stemming from the normal and anomalous 
loops, which we write here. The "normal" and "anoma- 
lous" terms of F xc for the electronic contributions, i 7 ^' 1 
and FSl> 2 , are as follows 



pei,i _ 



•i£(>-£)**>(i-£ 



tanh ( ^E k j tanh ( ^E k , 



F^ 2 = iV,(k,k')^^ 
A \ ' E k E k > 



kk' 



I anli ( ^E k j tanh (^E k , 



(50) 



(51) 



and the elcctron-phonon terms, with the normal and anoma- 
lous loops, F^c' 1 and FJ*' 2 i, respectively are given below 



pph,i 

xc 



\ J2 I dn a 2 F(Q) 

kk' J 



E k E k , 



I{E k ,E k ,,Q) 



pph,2 

xc 



1-^) I(E k ,-E k ,,f2) 

- Y [ df2 a 2 F(Q) 4^ 
2 ^ J E k E k , 



(52) 



x [I{E k ,E k ,,Q)-I{E k ,-E k ,,Q)]. (53) 
The function I(E k , E k >, Q) is defined as 

1 1 



I{E k ,E k ,,f2) 



-y 

fa 



(3 2 ^— ' iu>i - E k iu> 2 - E k i 
-2/2 



(ui - lu 2 ) 2 + S7 2 ' 
For the completeness, we give the definitions: 



(54) 



v{k, k') 



Aq 
9k,k+q 



a 2 F{f2) 



d 3 r d 3 r' ^(r)^k(r') 
1 



7T Vk'( r Vk'(r'), 
I-*- ^1 

J d 3 r ip k (r) V X q <y5k+ q (r), 
x 5(e k - e F ) <5(ek+q - £f), 



(55) 
(56) 

(57) 



where oj^q is the phonon frequency and N(ep) is the den- 
sity of states. 



Using the formulas (|50I54|) . one is ready to derive the 
exchange-correlation potential defined by Eq. (jl3|l . This 
derivation can be performed with the help of the chain 
rule as follows 



+ 



5F XC Sji 
SF XC S(a 



E 



SF XC 5\A 3 \ 2 



[5\A 3 \ 2 5 X * 



8(<t>j) 8 X * 



(58) 



Further evaluation of the above expression is given in de- 
tail in Refs. 5,12 . In this work, we give the final formula 
for A xc% i which involves the phonon and paramagnon spec- 
tral functions and can be implemented in a straightfor- 
ward way. We will give the details of implementation in 
Section EJ 

At this point, we arrived to the explicit expressions 
for F xc with the electronic and phononic parameters such 
as: the chemical potential fi, the density of states N(ef), 
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the single particle energies £k, and the Eliashberg function 
a 2 F(f2). Now, we are ready to introduce the spin fluctua- 
tions into the discussed formalism, and we will this in the 
following Section. 



5 Paramagnons in F xc [n, x] 

We will introduce the transverse spin-fluctuations to the 
total energy within the SCDFT. For the simplicity, we 
will assume the singlet pairing and the s-wave symmetry 
of the gap function. The extension to triplet supercon- 
ductors could be done following the work by Capelle et 
al. 16,20 . For the case of magnetic superconductors, one 
should take also into account a correction for the Zee- 
man effect, i.e. the spin gap. As for the paring potentials 
with the higher angular-momentum, one cannot average 
spherically the angular part of the interaction in the RPA 
formula for the paramagnon susceptibility, which formula 
will be used later in this Section. 

Here, we start with the Nambu Green's function for 
the superconductors with magnetic interactions included 
into the description. This matrix is now 4x4 dimentional 
and reads 



G(rr, rV) = -(TW\r, t) ® <P(r', /)), 



(59) 



with the 4-component field operators (the notation has 
been chosen according to Maki in Ref. ^Hl and x denotes 
the vector (r,r)) 



(60) 



The first-order selfenergy with the spin dependent in- 
teraction v 1 *" ', where \x and v denote the cartesian compo- 
nents of the spin orientations of two interacting electrons, 
is given by 



E{rr, rV) 



(rT,rV)a M G(rr,rVK, (61) 



'(rr,rV) = I ex (r)D^(r - T')I ex (r'). 



(62) 



The quantity I ex is the spin exchange interaction, and D^ u 
is the spin Green's function. The matrix is defined as 







(63) 



where a 1 ^ denotes a matrix transposed to the Pauli ma- 
trix (see Ref. HDD- 

For the transverse spin fluctuations, the a-matrix, given 
by formula l|63|l . involves the Pauli matrices er + and o~ 
defined as cr^ = \{cr x ± io- y ); explicitely 



1 







1 



Evaluation of the selfenergy with paramagnons, ac- 
cording to Eqs. H61I63|) . yields a very sparse 4x4-matrix 
which reads 



Z , (rr,r'r / ) = (rr, rV) 

/ G u (rr,r'r') -F iT (rr,rV)\ 


x 

y-F^rr.r'r'jOO G} T (rr, rV) j 



(64) 



where GL = — G^. 



Now, if we go back to the previous section and look 
again at the (l,l)-element of the (ZG)-matrix, we will 
remind to us that for the Coulomb and electron-phonon 
interactions, the total energy is proportional to the ex- 
pression l|46(l . For the magnetic interactions, however, for 
which the Nambu Green's function has been defined by 
Eq. I|59[) and the selfenergy has been given by Eq. I|64|) . 
the total energy is proportional to 



Gii Gtf — F\ tF-n — G 1 1 Gtt — Ft i F, 



(65) 



The above expression differs from relation (|46|) by signum 
in front of the anomalous Green's functions. This differ- 
ence will show up in the kernel Kij and the norm Zi of 
the gap equation such that, both the phonon and para- 
magnon spectral functions enter the kernel with different 
signum (originating from the anomalous loop of Green's 
functions) and the norm with the same signum (originat- 
ing from the normal loop). 

To proceed further with the evaluation of the xc-free 
energy, F xc , we write explicitely the spin- fluctuation Green's 
function, D^ v (t—t'), used in Eq. I|62|l . In the case of para- 
magnons, D^ v (t — t') is the transverse spin susceptibility, 
X H , defined as 

X + -(r-r',T-r') = (T S~{r, r)S+(r', r')>, (66) 

with the operators increasing and lowerring spin which are 
defined respectively as 



(67) 



S+(t,t) =$(r,r)^(r,T), 
S~(r,r) = ^{(r,r>0 T (r,T). 



For the conduction band, we can use a model of the 
homogeneous electron gas with the fluctuations treated on 
the level of the random phase approximation. The Fourier 
transform of the RPA-" dressed" paramagnon propagator 

x (q,O 



(69) 



with the Pauli susceptibility x° and the even Matsubara 
frequencies v n . 
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It is convenient to introduce the spectral representa- where the electronic part is defined by 



tion 



X +-(q,i/„) = - n^D°(n, 

J0 7T 



v n ) 3?m x + (q, 



K$ = toy, 
2tt 



D°(a,u n ) = 



-2Q 



(70) 
(71) 



log 



k'ikj \ ik% kj^j 2 ~\~ krp p 



(77) 
(78) 



and the momentum averaged paramagnon spectral func- 
tion 



P(Q) = N{e F ) 



X \I(q)f 



2kr 



dq 



<1 



2k p 



9m x + (?j ^) 



The Coulomb interaction uiij has been spherically aver- 
aged over the angular coordinates since, as we said before, 
we assumed the s-wave pairing. The electron correlations 
are taken into account by the Thomas-Fermi screening 
constant, kxF, and ki is an absolute value of the recipro- 
cal vector. 

The electron-phonon and -paramagnon interaction diago- 
nal part of the -matrix is given by 



(72) 



K 



ph+sf 



We assume that the interaction function, I(q), is the mo- 
mentum independent quantity I ex , which can be calcu- 
lated in a way given for instance in Ref. 

Therefore, for the systems with the electron-paramagnon 
interactions, the exchange-correlation free energy is given 

by 



_ 



kk' 
1 



1 - 



13 tanh(Pti/2)tanh(l3£j/2) 
x J dQ [a 2 F(Q) - P{Q)] 
x [m,Sj,n)-m (79) 
The diagonal part of the Mij -matrix is 

Zi = Zf + Zf +Sf , (80) 
where the purely electronic part is 



I(E k ,E k ,,Q) 
I(E k ,-Ev,Q) 



Zf 



(73) 



1 



^2wij[l-tanh{P^/2)} 



F 



s/,2 _ 



\ [ dQ p(q) 

kk' J 



E^Ev 



cosh' 2 (13 i k / 2) 

and the phononic and paramagnon part is 



[I{E k ,Eu,f})-I{E k ,-Eu,f2)], (74) 



where the function I(Ek,Ek>,Q) is defined by Eq. (|54*jl . 
The explicit formula for the paramagnon spectral func- 
tion, P(Q), within the RPA is given for instance in Refs. 



6 Gap equation with paramagnons and 
implementation details 

At this point, when we have completed the derivation of 
all components of the exchange-correlation free energy: 
the Coulomb part - Eqs. (|50I51|) . the phonon part - Eqs. 
H52I53I) . and the spin-fluctuation part - Eqs. (|73I74|) . we 
can write explicitely the gap equation given by Eqs. (|26t 



. (81) 



7 ph+sf 



-47T 



tanh(/3£i/2) (3 
x 2_, W2 sgn{u>2 



\ J dQ [a 2 F(Q) + P(Q)] 

. (82) 



r?ph + sf . r?ph + sf 



t,sym 



i,asym 



The My-matrix of the linearized equation (|26|l is the 
following function of the kernel Kij and the norm Zi 



1 Kjj[A = 0\ 
lJ 2 1 - Zi [A = 0] ' 

The nondiagonal part of the Mjj -matrix is given by 
K ■ - K el + K ph+sf 

J\ l3 — Ajj -|- J\ij , 



(75) 



(76) 



[w 2 2 + (& + Q) 2 ] 
+ +/„(&)] . 2 2( t~^I 21 2 - (83) 

+ (e t - ^) 2 ] 

Functions and are the Fermi-Dirac and Bose-Einstein 
distribution functions respectively. 

For the electronic part of the norm, i.e. Zf , we used 
the zero temperature approximation given in Refs. |14II10| . 
This approximation can be justified by the fact that the 
critical temperatures of simple metals, which we calculate 
in this work, are very low. The above simplification is done 
for sake of a numerical convenience since there are many 
singularities in the formula (|81(l . 

The subscripts "sym" and "asym" mean the symmetric 
and antisymmetric part of Zi with respect to the electron- 
phonon coupling elements <7k,k+q- The electron-paramagnon 
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interaction constant, /(q), has been also averaged in q 
leading to I ex . The antisymmetric part Zf^f^ is ommited 
in our calculations according to the reasons discussed in 
Refs. |5lll4) and in our previous work ^Hj. Therefore, we 
do not give the expression for Zf 1 ^^ in this work. 



7 Critical temperatures of simple metals 

In the following two subsections, we report the critical 
temperatures obtained by solving the SCDFT gap equa- 
tion with spin fluctuations included. We compare these 
results with the results without spin fluctuations and re- 
sults from the Eliashberg theory. First, we calculate pa- 
rameters of the gap equation for several simple metals: V, 
Mo, Ta, and Pd (fee and bec) at ambient pressure. At the 
end, we complete our previous results for Nb under pres- 
sure QUI reporting T c obtained within the SCDFT with 
the paramagnons included. 

The electronic structures, the densities of states (DOS) 
and the electron-phonon coupling constants and the phonon 
and magnon spectral functions for studied metals were 
calculated within the local density approximation (LDA). 
We used the pseudopotential plane wave codes pwscf 
|2T] and ESPRESSO |23 • The phonons and electron-phonon 
couplings were obtained from the density functional per- 
turbation theory (DFPT) 23 . Since the calculation of the 
spectral function a 2 F is very time consuming, we used 
the ultrasoft pseudopotentials (US PPs) [21] • The kinetic 
energy cut-offs for the wavefunctions and densities were 
45 Ry and 270 Ry respectively in order to reproduce well 
all features of the phonon dispersions especially for the low 
frequency phonons (see Ref. ^H])- The metallic broadening 
at the Fermi energy [22| was assumed at 0.03 Ry. We used 
the Monkhorst-Pack mesh [2H1 of (64,64,64)-points for the 
DOS calculations, (16,16,16)-points for the self-consistent 
calculation of the electron-phonon-coupling matrix ele- 
ments for the each phonon, the mesh of (8,8,8)-points 
to fit the phonon dispersions, and the fit from (16,16,16) 
into (64,64,64) mesh-points to perform the integrands with 
the double-delta function present in the definition of the 
electron-phonon coupling constant, X p , and the spectral 
function, a 2 F(uj). 

The spin-exchange interaction contants, I ex , for met- 
als at ambient pressure were taken from the work by Sin- 
galas et al. |23, an d further we used them for the calcu- 
lation of the spectral functions, P(lu), and the electron- 
paramagnon coupling constant, A sjF . For niobium under 
pressure, we used I ex and P{uj) calculated in our previous 
work [TUj . 

All electronic parameters and the phonon and magnon 
spectral functions were assumed to be the same for the 
normal and superconducting state. The accuracy of func- 
tions a 2 F(uu) and -P(w) is very important for an exact es- 
timation of the critical temperature. The electron-phonon 
spectral function, very time consumming for calculations, 
contains all the specific information about the studied sys- 
tem. In contrast to a 2 F(uj), the approximation which we 
used for the paramagnon spectral function, to avoid cal- 
culation of this quantity from the time-dependent density 



12 



C 10 
g 8 

a 

S? 6 
<D 
X> 




2 4 6 8 

frequency [THz] 

Fig. 1. The Eliashberg functions of V, Mo, Ta in bec structure 
and Pd in both bec and fee structures. 
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Fig. 2. The paramagnon spectral functions of V, Mo, Ta in 
bec structure and Pd in both bec and fee structures. 



functional theory, is insuficient. We made the assumption 
of the homogeneous electron gas for the spin susceptibility 
and the only spin-dependent quantity which we calculated 
specifically for a given metal was the exchange constant. 
The calculation of this constant, i.e. I ex , is very difficult 
and obtained results have a large error due to their very 
small values and necessity to calculate a response function 
to small magnetizations applied to the system. Therefore, 
as we will see below, the obtained critical temperatures 
are not always very close to the experimental ones. Fur- 
ther development should be directed into more accurate 
calculation of the spectral functions, especially P{uj). 



7.1 Transition metals at ambient pressure 

In TABLE ^ we report the critical temperatures and 
parameters which enter the gap equation calculated by 
means of the Eliashberg theory and the SCDFT for a 
few simple metals: vanadium, molibdenium and tanta- 
llum in bec lattice structure and palladium in fee and bec 
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Eliashberg SCDFT 

system N(e F ) \ ph X sf T c c p T c eps/ T c ep T c eps/ ~ T^ v "error" 



V bcc 


24.98 (26.14°) 


0.0218 6 


0.91 (1.19°) 


0.430 c 


Mo bcc 


8.81 (8.34°) 


0.0184 6 


0.47 (0.42°) 


0.024 c 


Ta bcc 


18.60 (18.38°) 


0.0162 6 


0.97 (0.86°) 


0.096 c 


Pd fee 


30.68 (34.14°) 


0.0230 6 


0.35 (0.35°) 


0.972 c 


Pd bcc 


16.60 (18.49°) 


0.0229 6 


0.68 (-) 


0.167 c 



0.212 


9.0 


5.9 


16.1 


7.4 


5.38 


38 


0.198 


0.8 


0.7 


1.5 


1.4 


0.92 


52 


0.209 


8.7 


8.1 


5.9 


4.6 


4.48 


3 


0.213 


0.01 












0.208 


1.3 


0.8 











° Values from Ref. [£B1 
b Values from Ref. $Z7\. 
c Calculated with I ex from Ref. |2"7) . 

Table 1. Various parameters such as: the crystal symmetry, density of states N(e_f) per Ry and per both spins, coupling 
constants lex [Ry/both spins], electron-phonon X ph , electron-paramagnon A s ^, and T c [K] calculated from the Eliashberg theory 
and the SCDFT with the Coulomb and phonon interactions only (ep) and with spin fluctuations (epsf), the experimental T^ xp 
(from Ref. W), and the "error" defined as (T c eps/ - T^ xp )/T^ xp [%] with T c eps/ calculated within the SCDFT. 



structures. Our calculated densities of states, N(ef), and 
electron-phonon coupling constants, X ph , are in a good 
agreement with previous calculations by Savrasov et al. 

The Eliashberg functions calculated within the DFPT 
are presented in FIG. ^ The Coulomb parameter, fj,*, 
was obtained from the Bennemen-Garland formula |29l 
I1U| . which employs the density of states. The spin ex- 
change constant, I ex , taken from Ref. has been used 
to obtain the paramagnon spectral function, P(u>), which 
we draw in FIG. |2] 

As for the critical temperatures, for tantallum, the 
SCDFT result is in a very small relative error, defined 
in TABLE ^ of 3% with respect to the experimental data 
pTT)] . While, the Eliashberg result with spin fluctuations 
included is in the error of 81%. For molibdenium, T c from 
the Eliashberg gap equation is smaller than the experi- 
mental one, even without the paramagnon effect. But the 
absolute error of all calculated temperatures for Mo is 
smaller than 1 K. Palladium in both structures fee and 
bcc is nonsuperconducting and the SCDFT reproduces 
well this result. In contrast to the SCDFT result, from 
the Eliashberg theory we obtained superconductivity for 
Pd in the bcc structure with a very small T c . 

Usually, the critical temperatures from the SCDFT are 
lower than temperatures from the Eliashberg theory. In 
some cases, however, the SCDFT temperatures are higher. 
This situation is for vanadium and molibdenium. Espe- 
cially for vanadium, T c from the SCDFT gap equation is 
about 2 K higher than the experimental data PU], even 
after inclusion of spin fluctuations. This fact may indicate 
that, either the spin exchange constant, J ex , was under- 
estimated, or a contribution of the asymmetric part of 
the phononic term in the SCDFT gap equation is quite 
large. As we know from results reported in Refs. [5TllO| . 
if we neglect the asymmetric part in the electron-phonon- 
coupling matrix elements by taking the a 2 F(ui) avaraged 
at the Fermi level, the critical temperatures are higher 
(see the discussion in SectionEJ. The last approximation, 



however, has to be done if we do not evaluate formulas 
with the .gk, k+q elements explicitcly. 

In general, the critical temperatures obtained from the 
SCDFT are in a good agreement with the measured tem- 
peratures |30| , and the effect of paramagnons improves the 
result considerable for many simple metals. 



7.2 Niobium under pressure 

In TABLE[3 we present critical temperatures and param- 
eters of the gap equation for niobium at eight pressures in 
the range from -17 GPa up to 80 GPa. The spin exchange 
constants, I ex , have been calculated from first principles 
in Ref. and the electron-phonon and electron-magnon 
spectral functions for Nb have been presented also in that 
work. 

Here, we complete our previous results by reporting 
the effect of paramagnons on T c calculated from the SCDFT. 
After the inclusion of spin fluctuations, the critical tem- 
peratures obtained from the SCDFT are closer to the ex- 
perimental T c 's for pressures in the range of 0-40 GPa, i.e. 
pressures between two anomalies measured by Struzhkin 
et al. 55 • The dependence of the measured critical tem- 
perature as a function of pressure is no longer reproduced 
by our calculations when we take into account param- 
agnons. At ambient pressure and for higher pressures, 
paramagnons seem to make the theoretical result worse. 
The above effect, could be explained by making the ob- 
servation that, in every case where the exchange constant 
I ex is large, the theoretical temperature underestimates 
the measured temperature, and vice versa, for the small- 
est I ex the critical temperature obtained from the SCDFT 
is the highest and the error is positive. 

Concluding this Subsection, the implementation of para- 
magnons to the SCDFT generally makes calculated criti- 
cal temperatures closer to the experimental ones. But our 
calculated exchange constants, I ex , are not sufficiently ac- 
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Eliashberg SCDFT 

N(e F ) I ex A ph A s/ T c ep T c eps/ T c ep T c eps/ ~ T^ p "error" 



-16.59 


22.82 


0.0211 


1.91 


0.28 


20.3 


16.7 


14.4 


6.2 






-9.45 


21.60 


0.0213 


1.60 


0.25 


19.5 


15.5 


13.2 


6.4 






-0.63 


20.24 


0.0217 


1.41 


0.22 


18.8 


14.7 


12.9 


7.2 


9.2 


-22 


9.98 


19.38 


0.0204 


1.65 


0.17 


19.6 


15.8 


13.4 


9.8 


10.0 


-2 


22.89 


18.32 


0.0189 


1.47 


0.13 


19.4 


16.0 


13.2 


11.3 


9.8 


15 


38.79 


17.10 


0.0228 


1.29 


0.16 


18.4 


14.1 


12.0 


10.1 


9.7 


4 


56.73 


15.42 


0.0292 


1.10 


0.23 


16.1 


10.7 


10.1 


8.4 


9.5 


-12 


78.37 


13.10 


0.0347 


0.86 


0.24 


13.7 


7.3 


8.2 


7.9 


8.8 


-10 



Table 2. Results for Nb; applied pressure p [GPa], density of states N(ef) per Ry and per both spins, spin exchange integral I ex 
[Ry/both spins] (from Ref. |10)L coupling constants: electron-phonon \ ph , electron-paramagnon A s ^, and T c [K] calculated from 
the Eliashberg theory (with /x*=0.21) and SCDFT with Coulomb and phonon interactions only (ep) and with spin fluctuations 
(epsf). The experimental T£ xp has been estimated from the picture given in Ref. |11|. Last column shows the "error" of the 
SCDFT calculations for T c eps/ defined in TABLE □ 



curate. This fact gives a direction for the future develop- 
ment. 



8 Summary 

In the present work, we included the transverse spin fluc- 
tuations to the density functional theory for supercon- 
ductors. The SCDFT is presented from its foundations, 
through the decoupling approximation, the gap equation 
and details of implementation. We assumed singlet and s- 
wave pairing potential; The extension to triplet supercon- 
ductors could be done following the work by Capelle et al. 
|lfiU2()| . The electron-phonon couplings and the electron- 
paramagnon couplings were averaged at the Fermi en- 
ergy, therefore the asymmetric part of the functional with 
respect to the electron-phonon matrix elements and to 
the spin-exchange interaction constants were ommited. 
Through the whole work, we kept the notation to be con- 
sistent with Parks [T7HT%| and Vonsovsky [T5|. 

Paramagnons and phonons in the superconducting state 
were assumed to be the same like in the normal state. The 
Eliashberg spectral function has been calculated within 
the density functional perturbation theory and it is fully 
material specific. Paramagnons, in contrast, have been ob- 
tained from the random phase approximation for the ho- 
mogeneous electron gas and only the spin exchange con- 
stants were calculated from the electronic structure. 

We reported the critical temperatures obtained from 
the SCDFT and the Eliashberg linearized gap equation 
with and without spin fluctuations for a few simple met- 
als: V, Mo, Ta, Pd at ambient pressure and Nb at several 
pressures up to 80 GPa. Some discrepancies between the 
temperatures calculated from the SCDFT and the mea- 
sured temperatures are due to the fact that it is quite 
difficult to obtain the accurate spin-exchange constants 
and/or to the fact that the spectral functions have been 
averaged at the Fermi level. Netherveless, the results show 



that inclusion of paramagnons improves the critical tem- 
peratures obtained from both methods, the SCDFT and 
the Eliashberg theory. The critical temperatures obtained 
from the parameter-free SCDFT are in most cases closer 
to the experimental data than the results obtained from 
the Eliashberg theory. 
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